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Using a simple relation between the virial expansion coefficients of the pressure and the entropy 
expansion coefficients in the case of the monomer-dimer model on infinite regular lattices, we have 
shown that, on hypercubic lattices of any dimension, the virial coefficients are positive through the 
20th order. We have observed that all virial coefficients so far known for this system are positive 
also on infinite regular lattices with different structure. We are thus led to conjecture that the virial 
expansion coefficients rrife are always positive. 

These considerations can be extended to the study of related bounds on finite graphs generalizing 
the infinite regular lattices, namely the finite grids and the regular biconnected graphs. The validity 
of the bounds A*^ln(i!A'^(i)) < 0 for k >2, where N(i) is the number of configurations of i dimers 
on the graph and A is the forward difference operator, is shown to correspond to the positivity of 
the virial coefficients. 

Our tests on many finite lattice graphs indicate that on large lattices these bounds are satisfied, 
giving support to the conjecture on the positivity of the virial coefficients. Moreover, in an exhaustive 
survey of some classes of regular biconnected graphs with a not too large number v of vertices, we 
observe only few violations of these bounds. We conjecture that the frequency of the violations 
vanishes as w —>■ oo. 

Using an inequality by Heilmann and Lieb, we find rigorous upper bounds on N{i) valid for arbi¬ 
trary graphs and for regular graphs. The similarity between this inequality and the one conjectured 
above suggests that one study the stricter inequality m*, > ^ for the virial coefficients, which is 
valid for all the known coefficients of the infinite regular lattice models. 

PACS numbers: 05.50.-|-q, 64.60.De, 75.10.Hk, 64.70.F-, 64.10.-|-h 
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I. INTRODUCTION 

The virial series-expansion coefficients of the pressure have been computed through relatively high orders for the 
monomer-dimer (MD) models on various infinite regular lattices. Presently, they are tabulated[T] through order 19 for 
the tetrahedral lattice and through order 7 for the hexagonal lattice[2]. Moreover, they are known for the triangular 
and the face-centered-cubic lattices through the orders 14 and 10 respectively [SH]. In the case of the linear lattice [5l6] 
and of the Bethe lattice [71H] . they are all known. We have recently computed [S] the expansions through the order 24 
for the hcc lattices in d = 3,4, 5,6, 7, and for the (hyper)-simple-cubic lattices, through order 24 in the case of the 
square, cubic and 4d lattices, through the orders 22 and 21 in dimensions d = 5 and 6 respectively, and finally [lO] in 
general dimensions d > 6 through the order 20. 

Long ago, Heilmann and Lieb |lllJ^ have studied the MD models also on finite graphs. They have shown that 
the zeroes of the matching generating polynomial M{z) — 'Y^N{i)z^ of a graph lie only on the real negative axis 
of the complex z plane. For the MD gas on a finite lattice, this implies the absence of any phase-transition in the 
thermodynamic limit. 

In this paper, we show that the first 20 coefficients of the virial expansion on hypercubic infinite regular lattices of 
any dimension are positive. With the knowledge that all the coefficients computed so far for the virial expansions of 
the MD models on all infinite regular lattices are positive, we are led to conjecture that all the virial coefficients are 
positive for the infinite regular lattice models. 

Using the definition of the graph dimer entropy [TJ], we argue that for a finite regular graph the bounds which 
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correspond to the positivity of the virial coefficients rrik for infinite regular lattices are 

A'=ln(d7V(i)) < 0 (1) 

with fc = 2, V and i = 0, v — k, where v is the matching number of G, i.e. the maximum number of pairwise 
disjoint edges of G. 

We have systematically tested the validity of Eq. Q within two classes of hnite graphs related to infinite regular 
lattices. The first class consists of the graphs induced by finite grids, the second consists of the biconnected regular 
simple graphs. The tests of Eq. Q for the latter class have been performed by exhaustively generating the biconnected 
regular graphs having a not too large number of vertices v, with the aid of the Nauty Drogram|15j via the Sage 
interface [T5] . 

The results of our survey are summarized in a set of schematic tables reported in the fourth Section. We have 
firstly restricted our attention to the class of finite lattices, being interested in finding sequences of finite lattices 
with no violations of the bounds in Eq. 0. This would indicate that all the virial coefficients are positive for the 
corresponding infinite regular lattices. Even though the bounds in Eq. Q have been derived only for regular graphs, 
it is however interesting to consider them also in the case of finite grids with open boundary conditions {be), (which 
clearly are not regular on the boundary). 

In summary, we found violations only for “narrow” lattices. This fact supports the conclusion that the bounds are 
not violated in the limit of inhnite lattices. 

We then turned to the more general class of the bipartite biconnected regular simple graphs. For the 3-regular 
graphs, we were able to test these bounds for u < 30 vertices. We have observed a few violations of the bounds for 

V > 18, but the frequency of these violations decreases regularly for v > 18. 

In the case of the bipartite biconnected 4-regular graphs, we have tested Eq. Q for u < 22. We have observed a 
single violation for v = 20. The frequency of the violations is even lower for v = 22. 

In the case of the bipartite biconnected r-regular graphs with r > 4, we could carry our tests only for v < 20, 
finding a single violation, that occurs in a 5-regular graph. 

When considering the non-bipartite biconnected 3-regular graphs, we observed few violations for v < 22, starting 
from V = 12; again the frequency of the violations decreases regularly after the hrst violation, for even values of 

V > 12. Notice that there are no such graphs with v odd. 

In the case of the non-bipartite biconnected 4-regular graphs, the frequency of the violations decreases regularly for 
even v after the first violation and a similar trend is observed for odd v. We could complete the tests only for v < 17. 
In this case one has to consider over 80 million graphs. 

Based on the results of our survey, we are led to conjecture that, for biconnected regular graphs, the bounds are 
violated with a frequency that vanishes as u —> oo. 

Let us now observe that the validity of the bounds of Eq. Q in the case k = 2 follows from the Heilmann-Lieb 
inequality Eq. (4) in |llj . Using this Heilmann-Lieb inequality, we shall derive rigorous upper bounds for N{i) of 
general graphs. In the case of regular graphs they are stricter than those in m in a region with small dimer density; 
in the case of general graphs, we obtain an upper bound for the matching generating polynomial which is stricter 
than the one found in [T8] . 

The (rigorous) Heilmann-Lieb inequality has the form 

A'=ln(d(n-z)LV(z)) < 0 (2) 

with n = [|] and k = 2. The similarity between this inequality and Eq. Q led us to investigate Eq. ([^ also for /c > 2 
and to conjecture for infinite regular lattices the stricter bound ruk > ^, which indeed is satished by all known virial 
coefficients. As in the previous case, we have investigated how well the bounds Eq. 0. with k > 2, are satisfied by 
finite grids and regular biconnected graphs. 

The tests of Eq. ([^, with k> 2, on lattice graphs give more violations than for Eq. Q. 

Therefore the examination of these lattices still gives some indication, though not as clear as in the case of Eq. Q , 
that the bounds Eq. 0 : with k >2, are satisfied by the corresponding infinite lattices and that the virial coefficients 
are positive. 

For the general graphs produced with the aid of Nauty, we find similar results as with the bounds of Eq. Q, 
although we observe more violations. 

We have also shown that, for any v, the bounds of Eqs. 00 are satisfied by the approximate distribution of 
bipartite regular random-graph obtained in [T5]. In the case of bipartite regular biconnected graphs, we argue that a 
conjecture on the entropy for these graphs, made in [Has], implies the conjecture that for u —>■ oo the frequency of 
the violations tends to zero. 

In all the tests performed on regular biconnected graphs (over more than 300 million graphs), Eq. 0 is valid for 
k = 3. We shall discuss also the corresponding upper bounds on N{i); the examples considered indicate that in the 
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case of bipartite regular biconnected graphs they are less strict than the Upper Matching conjecture (UMC) |19) . 
when the latter applies. 

Let us mention that while the UMC conjecture has not been yet proven, the lower asymptotic matching conjecture 
stated in m has been proven in [SD] and pT] . 

The paper is organized as follows. In the second Section, using a simple relation between the coefficients of the 
virial expansion of the pressure and the expansion coefficients of the dimer entropy, we prove that on hypercubic 
lattices the virial coefficients through order 20 are positive for generic d. In the third Section, we derive Eq. 0. 
obtain rigorous upper bounds for N(i) and discuss Eq. ([^. The fourth Section summarizes the results of the graph 
tests for a variety of lattices and graphs. In Appendix A we give the virial coefficients for d = 2,3 through order 
24. In Appendix B, the validity of Eq. Q, and of Eq. ([^ is proven for a few classes of graphs and for two average 
distributions. In Appendix C a formula for the bound on N(i) is deduced from the inequalities of Eq. 0. 


II. POSITIVITY OF VIRIAL COEFFICIENTS 


The combinatorial-statistical properties of a MD system on a r-regular lattice are usually described in the grand- 
canonical formalism of statistical mechanics, in which the pressure is defined as 


lim -ln5^(z) = P(z) = hz^. 

) —^OO 11 • ^ 


( 3 ) 


Here 'B.y{z) is the grand-canonical partition function for a u-site lattice and z is the activity. The dimer density is 
then 




2 = 1 


Setting p — 2/?, and solving Eq. 0 for z = z(p), we can express the pressure in terms of p 


P{p)=p/2 + Y, 


rrikP 


k=2 


This is the virial expansion. The entropy is defined by 

Kp) = -p(z)lnz -I- P(z) 

from Eqs.Q and Q one gets |10ll3] 

dX Inz 

dp 2 

Using the exDansion |10l2bj 

OO 

Hp) = R{p) + E 


A;=2 


where 


Rip) = 2 “ P^'^P ~ “ p)ln(l -p)-p) 


and r is the lattice coordination number, from Eq. Q one has 


In z = ln( 


r(I — pY 


- 2 ^ kakp'" ^ 


A:=2 


Substituting in Eq. ([^ Inz from Eq. (10) and A from Eqs.(|^|^ one obtains 

OO 

P = -ln(l - p) - ^ + ^(1 - k)akp’' 


( 4 ) 


( 5 ) 


( 6 ) 


( 7 ) 


( 8 ) 


( 9 ) 


( 10 ) 


fc=2 


( 11 ) 
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so that a simple relation is obtained between the coefficients of the virial expansion and the coefficients Ofc of the 
entropy expansion 


= {k - - ttk) ( 12 ) 

In the case of hypercubic lattices in any dimension d the coefficients Ofc have been computed in m through the 
order 20; in that reference in the cases d = 2,3,4 there are the values through order 24. In appendix A we list the 
corresponding virial coefficients for d = 2,3. 

On these lattices, we can express the coefficients of the pressure, of the virial and the entropy expansions as simple 
polynomials in the variable 1/d. Using the expressions for Uk with fc < 20 in Ref. [T(1] to examine the real roots of 
rrik, reported in Table [ij and observing that, for large d, the leading coefficient of mk is positive, it follows that rrik is 
positive for any integer d with d > 1. 

For example for k = 10 one has 


TOio 


1024d® - 35712d4 + 123240d3 - 118260d2 - 36990d + 67721 

10240d9 


with the three real root given in Table |Tj dividing the real axis in regions in which toiq has signs the 

positive integer values of d occur in the positive regions. 


TABLE I: Real roots of rrik for fc < 20 


k = 2 

0.25 




CO 

II 

-0.354 

0.354 



k — 4 

-0.859 




k — 5 

none 




II 

-0.239 




fc = 7 

-2.032 

0.848 



00 

II 

-1.796 

-0.557 

0.859 


it = 9 

1.044 

1.257 



fc = 10 

-0.655 

1.029 

1.313 


fe = 11 

-3.404 

0.998 



fe = 12 

-3.241 

-0.125 

0.997 


fc = 13 

1.000097 

1.725 



fc = 14 

0.9994 




fc = 15 

-4.657 

0.999997 



fc = 16 

-4.617 

1.000004 

1.801 


fc = 17 

1.000000085 

1.963 



00 

t-H 

II 

0.99999993 




fc = 19 

-5.852 

0.999999998 

2.005 

2.396 

fc = 20 

-5.879 

1.000000001 

1.993 



As remarked in the introduction, all the virial expansion coefficients so far computed for the lattice dimer models 
are positive. This leads us to conjecture that they are all positive on all infinite regular lattices. 

Let us remark that on any lattice with coordination number r one has 


^ 2k{k — 1) 


(13) 


for k less than the girth (i.e. the number of edges in a shortest cycle) of the lattice graph; in particular this is true to 
all orders in the case of the Bethe lattice m- 

Note that from the assumption that the virial coefficients are all positive, one gets an upper bound on the a^. 


1 

Ofc < 


k{k-l) 


(14) 











III. AN ARGUMENT FOR THE BOUNDS CONJECTURED IN EQ. @ AND EQ. 

In Ref. [14], we introduced the Newton series for the dimer entropy of a graph, in terms of the quantities 
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Here N{i) is the number of configurations of i dimers on the graph G with v vertices and N{i) given by 

v\ 




{v — 2i)!d2® 


(15) 


(16) 


is the number of configurations of i dimers on the complete graph on v = 2v vertices, where v is the matching number 
of G. If the graph G has a perfect matching v = v. 

For a graph that satisfies the “graph positivity” property introduced in Ref.[T3] one has 


A'=[d](f)>0 (17) 

where fc = 0,..., ^ and j = 0,..., r" — fc. The validity of these bounds for all allowed k and i is in fact equivalent to their 
validity for all allowed k and i = 0. 

If the “graph positivity” conjecture of Ref. |14] is true, then eq. 0 holds for almost all regular biconnected 
bipartite graphs (in a reasonable sense). More precisely, if this conjecture is true, then for each r the fraction of 
r-regular biconnected bipartite graphs with v vertices that satisfy eq. 0 tends to 1 as u —)■ oo. This positivity 
property is often satisfied also in non-regular bipartite graphs, while it is usually not satisfied by non-bipartite graphs. 

We have a similar situation here. Based upon the conjecture that all virial coefficients are positive for dimer models 
on infinite regular lattices, we state the following conjecture: the fraction of regular biconnected graphs with v vertices 
that satisfy 


A'' 


(1) < A'=ln 


{v — 2i)\v‘^^ 


(18) 


tends to 1 as u —^ oo, when k = 2,..., z/ and i = 0, ...,i> — k. The validity of these bounds for all allowed k and i is in fact 
equivalent to their validity for all allowed k and i = 0. Unfortunately, we do not know how to state mathematically 
how preponderantly this property holds for finite v. This will be seen “experimentally” in the following sections. We 


now trace the path which leads from positivity of the virial series coefficients to Eq. (181 in the case of finite lattices. 
With p « —, in the limit of large v one has (see Eq. (11) in [Uj l 


-d{i) akP^ 

V ^^ 


We eliminate ak using Eq. (12 ) and observe that 


1. (v — 

-In-^- 

V vl 


oo 


E 


k{k-l) 


(19) 


( 20 ) 


Thus we get 


1 ... 1 (u-2z)!t-2* “ 1 

-d{i) - -In^- ^ -^ > -— 

7? 7) 7?l I fC — 


k^2 


(k-1) 


nikP 


( 21 ) 


Using IA « (Uj, assuming that the are positive and that 1 for u —>■ oo, we are led to Eq. (181, in which 


we used v instead of v, motivated by the fact that the tests of Eq. (18) have fewer violations this way. The assumption 
that (( —> 1 is verified for biconnected 3-regular graphs and for regular bipartite graphs, since in these cases there 
exists a perfect matching (see |27j and references within). In fact all the biconnected regular graphs systematically 
examined in next section have v = n, where n = [|]. 
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Using Eqs.(15 16) we can rewrite Eq. (18) as 


A'^lniVU) < f-InU!) 

\ V — 1 

For k> 2, these bounds reduce to Eq. ([^. 

From the bound Eq. (4) in [TT], setting in that equation Zi = N{n — i) and M = n, with n = [|], we get 

(* + l)(^-*-l) _ a2 


( 22 ) 


A^lnA(i) < In- 


{i + 2){n — i) 


= —A^ln(d(n — i)!). 


(23) 


It follows that, for k = 2, the bounds of Eq. 0 are valid for any graph. 

Initially Eq. (18) would seem to hold only for the graphs, for example the periodic cubical graphs, whose limits 
are used to get the given lattice functions. But we will try to apply Eq. (18), or equivalently Eq. ([^, to regular 
biconnected graphs and to finite lattices. 


A. Rigorous upper bounds on the nnmber of matchings 


The bound Eq. (23) follows from the fact that all roots of the matching generating polynomial M{z) = Yli=o N{i)z’‘ 
are negativeso that the quantities Pv{i) = satisfy the Newton’s inequalities P{i + P)P{i — 1) < -P(*)^ j 

or equivalently 


I\hnP{i) < 0 


for i = 0,.., 1 / — 2. Using these inequalities and v <n, where n = [|], it is easy to see that also 


P{{) = 



(24) 


(25) 


satisfies them for i = 0,.., n — 2. 


Eq. (24) leads to the bounds (see [55] and Eq. (53) of Appendix C for the quantity g{i) = InP(i) with k = 2) 




V ^(*o) 

for * > * 0 - In the case of general graphs, we can apply Eq. 
is the number of edges of the graph. 


(26) 


with *0 = 0, using A(0) = 1 and iV(l) = E, where E 


In the following we will obtain a series of bounds based on Eq. (531 with given k and io; we will denote such bounds 


as “BXfc.io”, where B stands for Bound, X in G for a General graph, R for a Regular graph and B for Bipartite 
regular graph; k and ig are indices in Eq. (53). The case k = 2 corresponds to the Heilmann-Lieb inequality, holding 


for all graphs. The case fc = 3, discussed in the next subsection, is conjectural; to emphasize that we add the letter c 
to the name of the bound in the case k = 3. 

The choice of ig depends on the known N{i) for a category of graphs; for a general graph one has iV(0) = 1 and 
iV(l) = E. 

Thus we get an upper Bound for General graphs that we call BG2.0 to indicate that it follows from Eq. ([^ with 
k = 2,io = 0, 


Nii) < [E/nf 


From this a simple bound for the matching generating polynomial follows 


(27) 


\M{z)\ = \j2N{i)z^\<(l + ^)’ 

i=o \ ^ / 


(28) 


Using Pu{i) instead of P{i) one gets Eq. (28) with n replaced by i^; this stricter bound has been derived in [15] . 
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For E fixed, v < 2E. For the Hosoya index Z = M(l), introduced in theoretical chemistry to characterize the 
topological structure of large molecules, when v = 2E the bound Eq. }28| is saturated by the graph nK 2 . As u —)■ oo, 
for z = 1 the bound Eq. (28) tends to the bound Z < , first derived in [50] , 

In the following we will derive other bounds for P{i); one could write down analogous bounds using P^(i), when 
V < n. 

For a graph with E edges and vertices of degrees 5i, we have 



with 

V 

Y,S^=2E 


(29) 


(30) 


Define r = [^] and h = 2E — rv, then one has 0 < h < v. Defining 5i = r + ki, we get h = h and 


= vr"^ + 2r ki + > vr^ + (2r + 1) ki 

i—1 i—1 i—1 i=l 

so that 

V 

6i > vr^ + (2r + l)(2i? — rv) 

i=l 


(31) 


(32) 


It follows that using the bound Eq. 
for i > 0 


(26) with iQ 


1, we arrive at the bound BG2.1, slightly tighter than BG2.0 


N{i) < - 
n 


f 2Nm{2) Y ' 



where Nm{2) is the maximum value of N{2) 


Nra{2) = . 


E{E + 1) - - (2r + 1)(2E - rv) 


with 


r = 


-2E- 

- V - 


In the case of the BG2.1 bound, for the matching generating polynomial we obtain 


\M{z)\ < 1 + 


{n-l)E^ 

2nNra{2) 



2|z|iV„(2)\” 

{n-l)E) 


(33) 


(34) 


(35) 


(36) 


As an example, consider the logarithm of the Hosoya index Z for a general graph with u = 60 and E = 110. Then 
the upper bound BG2.0 for this quantity has the value 46.2 while the bound BG2.1 has the value 45.6. 

In [T7] it has been shown that for a regular graph with even u, the N{i) satisfy the following bounds for 0 < i < n; 
the bound 


N{i) < 2 -* 



(37) 


which is tight in the region of low dimer density and 


N{i) < 


(38) 











which is tight in the region of high dimer density. In |31| another bound tight in the region of high dimer density is 
given for i < n 


N{i) < exp ( -(pln(r) — _pln(p) — 2(1 — p)ln(l — p) — p + ) 

2 r — 1 


(39) 


where p = ^. We shall denote by BR, the bound for the regular graphs that is the minimum among these three 
bounds. For r constant and n large, the bound Eq. (pTl) improves Eq. (|37|) by a factor slightly larger than 


(1 -P) 


n-i+i 


( 1 _ |)2n-*+i 


for p not close to 0 or 1; for instance if i = this factor is 1.09”. Therefore a stricter bound (that we denote as 
BR2.0), is obtained combining Eq. (27) and Eq. (38). A slightly stricter bound (that we call BR2.1) is obtained 
combining Eq. (33)andEq. (38). 

In the case of regular bipartite graphs the following inequalities are known m for 0 < z < n 

N{i) < 

giving a tight bound in the region of low dimer density and 


(40) 


N(i) < 


(r!)- 


(41) 


giving a tight bound in the region of high dimer density. They lead to the bound Eq. (3.3) in [17] (that we denote 
as BB). The bound Eq. (27) is the same as Eq. (40) for these graphs. We can obtain a bound (denoted as BB2.2), 
which is stricter than BB, by using Eq. (26) with zq = 2, and observing that 

nr(rir — 2r + 1) 


and 


iV(3) = 


iV™(4) = 




+ 


N{2) = 


— 27z ( ) — nr{r — 1)^ — 2rz ( ) {nr — 3r + 2) 


2 2 
nr 


Using Eq. (26) with zq = 3, using fV(3) and the maximum value of fV(4), given by 

! 

-(l-2r) + 


-(19 — 60r + 52r^) + nr(- — 5r + 7r^ — 

24 4 2 ' 


rn(r — 1)^ 


(42) 


(43) 


(44) 


we get a bound slightly better than choosing zq = 2. For r fixed, v large, and zq = 3, the bound Eq. (26) for 7V(z) 
is smaller than the bound Eq. (40) by roughly a factor exp((3 — y)^^). For instance with z = f and n large, the 
bound is smaller by a constant factor exp(|(l — y)). We denote by BB2.3 the bound Eq. (26) for Zq = 3 combined 
with the bound Eq. (41). 


B. Conjecture of a stricter bound on virial coefficients 

We now set p = p{i) = 2ijv and assume that v, n — i are large. From Eq. ( |21|16|15[ ) it follows that for k > 2 

nijp^ 


A'=(-ln(iV(z)z!) « 


i=2 


j - 1 


(45) 


from which it follows that 
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All the known virial coefficients in the infinite regular lattice models mentioned in the introduction satisfy the bound 




(47) 


for fc > 2. We have also checked that 
satisfy the bound Eq. 
these cases. We are therefore lead to conjecture that Eq. 


(47) through order 20 for any d, so that from Eq. 

(47) holds on a) 


in d dimensional hypercubic infinite lattices the virial expansion coefficients 

(46) it follows that Eq. (§ is satisfied in 
1 infinite regular lattices. This indicates 
for p = 1. As already indicated in [10] the entropy 
1. However from Eq. ([^ if follows that ^ is 


that there is a singularity in the virial expansion of the pressure, 
is analytic in the interval 0 < p < 1 and it is bounded above at p = 
log-singular for large z (equivalently for p —>■ 1). As to the pressure, from Eq. and the fact that A(p) is bounded 
above at p = 1 and p(z) —>■ 1 as z —?► oo, it follows that the pressure is log-singular for large z. 

Eq. (U) for fc = 2 is satisfied for any graph, due to the Heilmann-Lieb inequality Eq. ( |2^. We are therefore led to 
the conjecture that in all infinite regular lattice models the virial coefficients satisfy Eq. (47), and to the conjecture 
that the frequency of the violations to the bounds Eq. ([^ for large regular biconnected graphs, tends to zero as 

V —t 00. 

In [19] , |25j it has been conjectured that in the limit of infinitely large random bipartite regular graphs, the entropy 
has almost surely the coefficients ak given in Eq. (13); the same is true for the average distribution Eq. (49) introduced 
in [19]. In the appendix we show that this average distribution satisfies the bounds Eq. ([^, so that this conjecture 
suggests that the bipartite regular graphs almost surely satisfy Eq. 

In the tests discussed in the next section, no violation of the bound Eq. ([^ with fc = 3 are observed both for regular 
biconnected graphs and for not necessarily regular lattice graphs. Notice however that for fc = 3, in the case of more 
general non-regular graphs there are violations, e.g. in the case of the graph with perfect matching shown in Eigure 


m 


It seems therefore interesting to investigate the bound Eq. ([^ for fc = 3 in the case of biconnected regular graphs. 
Using Eq. (53) for fc = 3 we get 


P(z) < P{io) 


Pjio + 1 ) 
Pi^o) 


fP{io + 2)P{io)\ 

\ Pi^o + ir ) 




(48) 


for z > iQ. 

In the case of regular graphs, we can use Eq. 
conjectured bound BRS.Oc is obtained combining this bound with Eq. (38). 


(48) with zo = 0, taking into account Eq. (42) for N{2). The 


As an example of a regular graph, consider the Buckminster fullerene Cgo; it has u = 60 and E = 90; the first 
violation of Eq. ([^ is for fc = 20. The logarithm of the Hosoya index Z of this graph has the value ln(Z) = 34.89. 
The bound BR has the value 46.49, the bound BR2.0 yields 41.50, the bound BR2.1 is 41.02, while BR3.0c is 36.58. 

In the case of regular bipartite graphs, we can use Eq. (48) with zq = 2, using N{2),N{3) and replacing iV(4) by 
A^m(4), as in subsection BIA. We thus conjecture a bound (railed BB3.2c) obtained combining this bound with Eq. 


(41). 


When r divides zz, the bound BB3.2c is weaker than the Upper Matching conjecture (UMC)[TS], according to which 
the number of z-matchings of a r-regular bipartite graph with n = qr is bounded above by the number of z-matchings 
of qKr^r, for q a positive integer. 

In Table jll] we give InZ, where Z is the Hosoya index, and upper bound estimates for two graphs, the periodic 
12 X 12 square grid and the 6-cube; for the latter we took the matching polynomial from [^ which used [23124] . 


TABLE II: \nZ, where Z is the Hosoya index and upper bounds estimates for two bipartite regular graphs. The first graph is 
the periodic grid 12 x 12, with V = 144 and degree 4; the second one is the 6-cube, with U = 64 and degree 6. The latter two 
upper bounds BBS.2c and UMC are conjectural. The UMC conjecture applies only to graphs in which U/2 is a multiple of the 
degree, so that for the 6-cube no entry appears in the Table. 


graph 

In^ 

BB 

BB2.3 

BB3.2c 

UMC 

periodic grid 12 x 12 

6-cube 

95.44 

50.32 

115.24 

59.59 

113.69 

58.49 

100.11 

52.87 

96.16 


The d-cubes through d = 6 satisfy the bounds Eq. ([^ and Eq. ([^. 
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FIG. 1: A graph violating the bound Eq. § for k — S. 


IV. TESTS ON THE UPPER BOUNDS 

We shall now review our graph tests. All graphs we considered are simple biconnected graphs. The tests of Eqs.Q 
and Eqs.([^ are done for all allowed values of i (one could check only the case t = 0 to see if there are violations for 
some k, but the value of k for which there is a violation is generally higher if one considers only the case f = 0). 

We made tests on lattice graphs and on regular graphs; to generate systematically the latter we have used the geng 
program in the Nauty package [IS], via the Sage interface|16j. The matching generating polynomials are computed 
with the aid of the algorithm described in [33]. To perform our computations, we have used an ordinary desktop 
personal computer based on a processor Intel i7 860 with a RAM of 8 GB. 
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A. Tests of the bounds Eq. (Q on finite lattice graphs, i) Periodic boundary conditions 

Let us first discuss the tests on finite lattices with periodic be. In the case of rectangular grids of size x Ly with 


Lx > Ly, we considered the cases listed in Tab. Ill 


TABLE III: Violations of the bounds Eq. (|^ in the case of rectangular grids of size Lx x Ly , with Lx > Ly and periodic be. 


Lx 

Ly 

Violations of the bounds Eq. ( 

2 

1 

Lx < 2000 

Ly = 3 

for Lx >9 

Lx < 1100 

II 

for all Lx > 409 

Lx < 500 

Ly ~ 5 

for Lx > 186 

Lx < 200 

Ly = 6 

none 

Lx < 70 

Ly=7 

none 

Lx < 50 

Ly = 8 

none 

Lx < 20 

Ly = 9 

none 

Lx < 15 

Ly = 10 

none 

Lx = 11 

Ly = 11 

none 

Lx = 12 

Ly = 12 

none 


Therefore we have found violations only for Ly < 5, when Lx > 3Ly. The minimum k for which there are violations 
is 12 for the cases with Ly = 3, it is larger than 100 in the other cases. 

In the case of triangular grids of size Lx x Ly with Lx > Ly, we have considered the cases listed in Tab. 


TABLE IV: Violations of the bounds Eq. 0 in the case of triangular grids of size Lx x Ly , with Lx > Ly and periodic be. 
These grids are obtained by adding a SW-NE diagonal to rectangular grids Lx x Ly. 


Lx 

Ly 

Violations of the bounds Eq. ( 

2l 

Lx < 1000 

Lx < 1000 

Lx < 400 

Lx < 200 

Lx < 50 
Lx < 17 

Ly = 3 

L, =4 

Ly = 5 
Ly = 6 

Ly = 7 

Ly = 8 

for all Lx > 16 

none 

for Lx > 66 

none 

none 

none 


IV 


In the case of hexagonal lattices of size Lx x Ly with periodic be (in the brick-wall representation) with Lx, Ly even, 
we have considered the cases listed in Tab. El 

TABLE V: Violations of the bounds Eq. 0 in the case of hexagonal grids of size Lx x Ly , with Lx > Ly and periodic he. 


Lx 

Ly 

Violations of the bounds Eq. ( 

2 

1 

II 

4 < Ly < 100 

all cases 

Lx = 6 

4 < Ly < 100 

for Ly > 26 

8< Lx<U 

4 < Ly < 14 

none 

6<Lx< 100 

Ly = 4 

none 

6<Lx< 100 

Ly = 6 

none 


The minimum value of k for which we found a violation is 7. In particular there are no violations for Lx close to 


Ly and Lx > 4. 


B. Tests of the bounds Eq. (Q on finite lattice graphs, ii) Open boundary conditions 

Let us now turn to finite lattices with open boundary conditions. 

In the case of rectangular grids of size Lx x Ly, with open be, we have examined the cases 2 < Ly < Lx < 19, 
finding no violations to Eq. 0. 
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In the case of triangular grids of size x Ly, obtained by adding a SW-NE diagonal in a rectangular grids x Ly 
with open be, for 2 < Ly < < 18 there are violations only for Ly x 3 with L^ > 10. The minimum value for which 

Eq. 0 is violated is k = 16. 

We have also examined three-dimensional slabs of size Lx x Ly x L^ with open be. The results are summarized in 

Tab. ED 


TABLE VI: Violations of the bounds Eq. 0 in the case of three-dimensional slabs of size L^ x Ly x with open he. 


Lx 

Ly 

L. 

Violations of the bounds Eq. ( 

2 

1 

2<Lx< 1000 

Ly = 3 

L, = 2 

none 

3< Lx< 500 

11 

CO 

II 

CO 

none 

A< Lx< 600 

II 

II 

to 

none 

A< Lx< 200 

II 

L, = 3 

none 

4 < < 40 

Ly — 4 

II 

N 

none 


For all the graphs examined in this section, Eq. 0 is satisfied for k < A. It would be interesting to know whether 
these bounds, Eq. 0 for k < A, are always satisfied for regular biconnected graphs. 

C. Tests of the bounds Eq. (0 on bipartite graphs 


We have tested the validity of the upper bounds Eq. (181 (equivalently Eq. 0 ) for regular bipartite biconnected 
graphs (RBB), by enumerating and studying exhaustively a large class of graphs. For u < 18 vertices, we observe a 
single violation in the case of a 3-regular graph with v = 18. For 3-regular graphs with 18 < u < 30, the frequency of 
the violations decreases with increasing v (v is even, since there are no graphs with v odd in this class). 

It is interesting to observe that, in the cases considered in Table E^D the average order of the automorphism groups 
of the positivity-violating graphs is a few times larger than the average order of the automorphism groups of all the 
RBB graphs with the same vertex degree. The same is true in the following tests of this section. 


TABLE VIL For the RBB graphs with a given number 18 < n < 30 of vertices of degree 3, we have listed the number of graphs 
in this class, the number of violations of the upper bounds, the average order ng of the automorphism groups of the graphs, 
the average order ngv of this group for the graphs violating the bounds Eq. 0 - k is the minimum value for which Eq. 0 is 
violated. 


V 

number of graphs 

violations 

ng 

ngv 

k 

18 

149 

1 

15.1 

64. 

9 

20 

703 

3 

8.7 

91. 

10 

22 

4132 

13 

4.5 

40. 

9 

24 

29579 

38 

3.3 

32. 

7 

26 

245627 

253 

2.3 

22. 

7 

28 

2291589 

1392 

1.9 

20. 

6 

30 

23466857 

8008 

1.7 

16. 

7 


For the RBB graphs with vertices of degree 4 there is one violation for v = 20 among 62611 graphs and fc = 10 is 
the minimum value for which Eq. 0 is violated. The order of the automorphism group of the violating graph is 256, 
while the average order is 3.1. For v = 22, there are 5 violations among 2806490 graphs and the minimum value for 
which Eq. 0 is violated is k = 11; the average of the orders of the automorphism groups of the violating graphs is 
3721.6, while the average order of all graphs is 1.5. 

For the RBB graphs with vertices of degree 5, there is one violation for v = 20 out of 304495 graphs, with k = 9; 
the order of the automorphism group of this graph is 1327104, while the average order is 7.1. 

For the RBB graphs with vertices of degree larger than 4, we had to restrict to graphs with w < 20 vertices and 
observed no violations for degree larger than 5. 
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D. Tests of the bounds Eq. ([^ on non-bipartite graphs 

In the case of the biconnected 3-regular non-bipartite graphs with v < 22, the first violation occurs for v = 12. For 
u > 12 the frequency of the violations decreases with v, as shown in Table |Vni| 


TABLE VIII: For the biconnected 3-regular non-bipartite graphs with a given number v of vertices, we have listed the number 
of graphs in this class, the number of violations of the upper bounds, the average order ng of the automorphism groups of the 
graphs, the average order ngv of this group for the graphs violating the bounds Eq. 0- The cases with v odd are not listed, 
since there are no such graphs, k is the minimum value for which Eq. 0 is violated. 


V 

number of graphs 

violations 

ng 

ngv 

k 

12 

76 

1 

7.4 

16. 

6 

14 

467 

6 

4.4 

8.7 

6 

16 

3836 

44 

3.1 

8.3 

5 

18 

39717 

257 

2.2 

7.5 

5 

20 

497115 

2856 

1.7 

5.6 

5 

22 

7183495 

29597 

1.5 

4. 

5 


In the case of the biconnected 4-regular non-bipartite graphs with v < 17, the first violations occur for v = 12, as 
shown in Table Hxl 


TABLE IX: For the biconnected 4-regular non-bipartite graphs with a given number v of vertices , we have listed the number 
of graphs in this class, the number of violations of the upper bounds, the average order ng of the automorphism groups of the 
graphs, the average order ngv of this group for the graphs violating the bounds Eq. 0 . k is the minimum value for which Eq. 
Q is violated. 


V 

number of graphs 

violations 

ng 

ngv 

k 

12 

1538 

2 

3.4 

40 

6 

13 

10768 

0 




14 

88112 

12 

1.6 

17 

6 

15 

805281 

30 

1.3 

14. 

7 

16 

8036122 

454 

1.2 

14. 

6 

17 

86214189 

295 

1.2 

10. 

6 


Unlike in the bipartite case and in the case of 3-regular non-bipartite graphs examined above, there exist graphs 
with odd V. 

The frequency of the violations occurring among the graphs with even v decreases regularly as v increases. The 
same is true for the graphs with odd v. 


E. Tests of the bounds Eq. ([^ on finite lattice graphs, i) Periodic boundary conditions 

Let us first discuss the tests on finite lattices with periodic be. 

In the case of rectangular grids of size x Ly with > Ly, we have considered the cases listed in Tab. 
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TABLE X: Violations of the bounds Eq. § in the case of rectangular grids of size x Ly , with Lx > Ly and periodic be. 


Lx 

Ly 

Violations of the bounds Eq. ( 

2 

1 

Lx < 2500 

II 

CO 

for all Lx > 5 

Lx < 1000 

II 

for most Lx >8 

Lx < 700 

Ly — 5 

for most Lx > 9 

Lx < 200 

II 

for most Lx > 12 

Lx < 70 

Ly — 7 

none 

Lx < 50 

II 

00 

for most Lx > 14 

Lx < 20 

Ly ~ 9 

none 

Lx < 15 

Ly = 10 

none 

Lx = 11 

Ly = 11 

none 

Lx = 12 

Ly = 12 

none 


The minimum value of k for which violations are observed, is 4. 

In the case of the triangular grids of size x Ly with > Ly and periodic 6c, we considered the cases listed in 

TabKH 


TABLE XL Violations of the bounds Eq. ([^ in the case of triangular grids of size Lx x Ly , with Lx > Ly and periodic be. 


Lx 

Ly 

Violations of the bounds Eq. (2 1 

Lx < 2500 

Lx < 1000 

Lx < 500 

Lx < 150 

Lx < 50 

Lx < 15 

Ly = 3 
Ly = 1 
Ly = 5 
Ly = 6 
Ly =7 
Ly = 8 

for most Lx > 11 

for all Lx > 88 

for most Lx > 29 

none 

none 

none 


The minimum value of k for which there are violations is 4. 

In the case of the hexagonal grids (in the brick-wall representation) of size L^ x Ly, with L^ > Ly and periodic be 
we considered the cases listed in Tab IXIIl 


TABLE XIL Hexagonal grids (in the brick-wall representation) of size Lx x Ly, with Lx > Ly and periodic be tested . 


Lx 

Ly 

violations of the bound Eq. ( 

21 

Lx = 4 

Lx = 6 

8< Lx<l4 

8<Lx< 100 

8< Lx <30 

4 < Ly < 100 
4 < Ly < 100 
4 < Ly < 14 

Ly=A 

Ly = 6 

Ly = 6 and Ly >= 10 

Ly >0 

none for Lx — 8 < Ly < Lx 

Ly > 14 

Lx > 18 


Summarizing this table, we found no violations in the band max{4,, Lx — 8) < Ly < Lx', almost all the other cases 
violate the bounds. The violations of the bounds Eq. ([^ occur for k > 4. 


F. Tests of the bounds Eq. ([^ on finite lattice graphs, ii) Open boundary conditions 

Let us now turn to finite lattices with open be. 

In the case of the rectangular grids of size Lx x Ly with 2 < Ly < Lx < 19 and open be, there are many violations, 
but they do not occur for Lx = Ly with Lx ^ 5, or for Ly > 9. The minimum value of k for which Eq. ([^ is violated 
is fc = 4. 

In the case of the triangular grids Lx x Ly with 2 < Ly < Lx < 18 and open be, there are no violations for Lx = Ly 
with Lx ^ 5, or for Ly > 8. The minimum k for which there are violations to Eq. (j^ is A: = 6. 

In the case of the 3-d grids of size Lx x Ly x L^ with open be, we considered the cases listed in Tab. 


XIII 
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TABLE XIII: Violations of the bounds Eq. § in the case of three-dimensional slabs of size Lx x Ly x Lz with open be. 


Lx 

Ly 

Lz 

Violations of the bounds Eq. ( 

b 

Lx < 1000 

II 

CO 

Lz=2 

for Lx >3 

Lx < 500 

CO 

II 

Lz=3 

for Lx = 5,9,11,13,15,17,21 

Lx < 600 

Ly ~ A 

Lz=2 

for Lx > 11 

Lx < 200 

Ly = A 

Lz=3 

for Lx > 16 

Lx < 40 

II 

Lz=A 

for Lx > 20 


In particular, no violations are observed for = Ly = L^. 

Summarizing, in the tests of Eq. on lattice graphs, violations are observed more often than for Eq. 0; for 
bidimensional lattices with ~ Ly > h there are no violations, except in the case of the hexagonal lattice, in which 
no violations are present on a strip with > Ly. This restriction might be due to the fact that the hexagonal lattice 
is sensitive to the boundary conditions[34j. Therefore these tests give some indication, although not as sharply as in 
the case of Eq. that in the limit of infinite lattices the bounds Eq. 0 and the corresponding bounds Eq. ( |47[ ) 
are satisfied, and virial positivity holds. 


G. Tests of the bounds Eq. (0 on bipartite graphs 


Let us now discuss the systematic tests on regular graphs. 

In the case of the RBB graphs with vertices of degree 3 and v < 30, the first violation occurs for v = 14; for larger 
values of v the frequency of violations decreases irregularly (e.g. it increases a.t v = 20), as shown in Table XIV 


TABLE XIV: For the RBB graphs with a given number 14 < v < 30 of vertices of degree 3, we have listed the number of 
graphs in this class, the number of violations of the upper bounds Eq. 0 , the average order ng of the automorphism groups of 
the graphs, the average order ngv of this group for the graphs violating the bounds Eq. 0; A: is the minimum value for which 
these bounds are violated. 


V 

number of graphs 

violations 

ng 

ngv 

k 

14 

13 

1 

44. 

28. 

7 

16 

38 

2 

19. 

48. 

6 

18 

149 

5 

15. 

84. 

5 

20 

703 

33 

8.7 

29. 

5 

22 

4132 

106 

4.5 

20. 

4 

24 

29579 

619 

3.3 

16. 

4 

26 

245627 

3415 

2.3 

10. 

4 

28 

2291589 

22913 

1.9 

7.9 

4 

30 

23466857 

163789 

1.7 

6.2 

4 


In the case of the RBB graphs with vertices of degree 4 and v < 22, the first violation occurs for v = 16. Eor larger 
values of v the frequency of the violations decreases, as shown in the Table |XVj 

In the case of RBB graphs with vertices of degree 5 and v < 20, the first 3 violations occur for v = 20 out of 
304495 graphs, with average order 449280 of the automorphism groups of the violating graphs, to be compared with 
an average order of 7.1 for all the graphs; the minimum k for which there are violations is A: = 6. 

We have checked that for RBB graph with v < 20, there are no violations for v > 5. 


H. Tests of the bounds Eq. 0) on non-bipartite graphs 


In the case of biconnected 3-regular non-bipartite graphs with v < 20, the first violation occurs for v = 10. For 
larger values of v the frequency of violations decreases irregularly as shown in Table |XVI| 

In the case of biconnected 4-regular non-bipartite graphs with v < 17, a first violation is met for r; = 11. For larger 
V, the frequency of the violations decreases regularly for v even, while it decreases irregularly for v odd (it increases 
for z; = 15), see Table 


In addition to the systematic examination of 3-regular biconnected graphs up to ?; = 22, we have tested the k = 3 
bound for v = 30. Using the NetworkX [35j random regular graph generator, we have thus examined over 100 millions 
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TABLE XV: For the RBB graphs with 16 < v < 22 vertices of degree 4, we have listed the number of graphs in this class, 
the number of violations of the upper bounds Eq. it- the average order ng of the automorphism groups of the graphs, the 
average order ngv of this group for the graphs violating the bounds Eq. ([^; fc is the minimum value for which these bounds 
are violated. 


V 

number of graphs 

violations 

ng 

ngv 

k 

16 

129 

1 

112. 

5184. 

8 

18 

1980 

1 

8.7 

576. 

7 

20 

62611 

18 

3.1 

1901. 

6 

22 

2806490 

115 

1.5 

487. 

6 


TABLE XVI: For the biconnected non-bipartite graphs with a given number 10 < w < 22 of vertices of degree 3, we have 
listed the number of graphs in this class, the number of violations of the upper bounds Eq. it- the average order ng of the 
automorphism groups of the graphs, the average order ngv of this group for the graphs violating the bounds Eq. (|^; A: is the 
minimum value for which these bounds are violated. 


V 

number of graphs 

violations 

ng 

ngv 

k 

10 

16 

1 

15. 

16. 

5 

12 

76 

2 

7.4 

12. 

5 

14 

467 

11 

4.4 

11. 

4 

16 

3836 

102 

3.1 

9. 

4 

18 

39717 

741 

2.2 

6.2 

4 

20 

497115 

7324 

1.7 

4.8 

4 

22 

7183495 

78426 

1.5 

3.5 

4 


3-regular graphs with v = 30, and have checked that Eq. ([^ is satisfied for fc = 3. From Table XVI by a simple 
extrapolation we estimate that for u = 30 there are roughly 8 x 10^^ graphs, so the fraction of non-inequivalent 
produced random graphs is expected to be less than 10“^. 

Analogously, we have examined 100 millions r andom 4-regular biconnected graphs with 21 vertices. We have checked 
that Eq. ([^ is satisfied for fc = 3. From Table XVII we estimate that for v = 21 there are roughly 2 x 10^^ graphs, 
so the fraction of equivalent random graphs produced by the generator, is expected to be less than 10“^. 

In all the tests performed on regular graphs (over more than 300 million graphs) no violations of Eq. (H for A: = 3 
are observed. 

As a further comment, we observe that there are some similarities between the graph positivity property J3] and 
Eq. 0: in both cases the violating graphs have an average order of the automorphism groups which is several times 
larger than the average over all the graphs with the same number of vertices v. This property is observed also in the 
case of Eq. 0, but to a lesser extent. 


V. CONCLUSIONS 


We have observed that the coefficients of the virial expansions of the MD model computed up to now satisfy Eq. 


(47) on all infinite regular lattices. In particular the first 20 coefficients of the virial expansion satisfy Eq. (47) for the 


hypercubic infinite lattices of any dimension d. This led us to conjecture that the virial coefficients are all positive for 
any infinite regular lattice model, and to the stricter eonjecture that they all satisfy Eq. (47). 


Using a simple relation between the virial coefficients and the coefficients of the series for the dimer entropy, the 
eonjecture on the positivity of the virial coefficients led us to test the validity of the bounds in Eq. Q for the finite 
lattice graphs and also for the finite graphs which somehow generalize them, namely the biconnected regular graphs. 
For k = 2 these bounds follow from the inequality Eq. (23) proved by Heilmann and Lieb. 


We have shown that this inequality leads to rigorous upper bounds for the number of matchings N(i) improving 
those known up to now for regular graphs [32] in the region of low dimer density, and for general graphs. In the latter 
case, we derive an upper bound for the matching matching polynomial improving the one in m- 

The tests on lattice graphs support the validity of the virial positivity conjecture. 

Our tests on a large class of regular graphs and lattice graphs also support the conjecture that the frequency of 
violations of Eq. Q tends to zero as u —)■ oo. 

We have proposed the more general bounds Eq. 0 corresponding to the stricter conjecture that Eq. ( |47| is valid 
for infinite regular lattices . 
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TABLE XVII: For the biconnected non-bipartite graphs with 10 < n < 17 vertices of degree 4, we have listed the number of 
graphs in this class, the number of violations of the upper bounds Eq. § , the average order ng of the automorphism groups of 
the graphs, the average order ngv of this group for the graphs violating the bounds Eq. ([^; A: is the minimum value for which 
these bounds are violated. 


V 

number of graphs 

violations 

ng 

ngv 

k 

11 

264 

1 

5.2 

12. 

5 

12 

1538 

3 

3.4 

32. 

5 

13 

10768 

37 

2.1 

16. 

5 

14 

88112 

34 

1.6 

31. 

5 

15 

805281 

3086 

1.3 

5.0 

4 

16 

8036122 

1121 

1.2 

14. 

5 

17 

86214189 

197431 

1.2 

2.5 

4 


Our tests of the bounds Eq. ([^ on the regular graphs give results similar to those obtained for the bounds Eq. Q , 
although with more violations. 

The tests of Eq. (|^ on lattice graphs show many more violations than for the bounds Eq. 0 ; there are however 
bands with « Ly in the 2d case and with ~ Ly for the cubic slabs, for which there are no violations. 

Extrapolating the behavior in these regions to the limit of large lattices we get some indication that there are no 
violations of bounds Eq. ([^ and the inequality Eq. ( |47[ ) in the infinite lattice limit, but this indication is not as 
strong as in the case of the bounds Eq. 0 and virial positivity. 

In all tests (carried out over more than 300 million regular graphs), we found no violations of the bounds Eq. 0 for 
k = 3. It would be interesting to know whether the bounds Eq. 0 for /c = 3 are satisfied by all regular biconnected 
graphs. 

We have discussed upper bounds that could be proved if Eq. 0 held for fc = 3: they are tighter than those we 
proved based on the truth of the k = 2 case. 


VI. APPENDIX A: VIRIAL COEFFICIENTS IN d = 2, 3 


The virial coefficients ruk in d = 2,3,4 can be obtained from the entropy coefficients computed in m through 
order 24; they satisfy the bound Eq.(47|. In Table ?? we give the virial coefficients in d = 2, 3. 
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TABLE XVIII: Virial coefficients mk in d = 2, 3 up to order 24. 


k 

d = 2 

d = 3 

2 

7/16 

11/24 

3 

31/96 

71/216 

4 

121/512 

419/1728 

5 

471/2560 

31/160 

6 

1867/12288 

15031/93312 

7 

7435/57344 

89951/653184 

8 

29477/262144 

539963/4478976 

9 

116383/1179648 

3244127/30233088 

10 

459517/5242880 

19482611/201553920 

11 

1821051/23068672 

116960471/1330255872 

12 

7255915/100663296 

702028151/8707129344 

13 

29063919/436207616 

1404544085/18865446912 

14 

16697149/268435456 

8431212005/121899810816 

15 

157001097/2684354560 

151861682911/2350924922880 

16 

1898046421/34359738368 

911893249819/15045919506432 

17 

7634823999/146028888064 

5475822286271/95917736853504 

18 

30619146937/618475290624 

32879301057475/609359740010496 

19 

122399296903/2611340115968 

197399821995527/3859278353399808 

20 

488028559661/10995116277760 

395010395194633/8124796533473280 

21 

1943383170991/46179488366592 

2371215487550117/51186218160881664 

22 

7740357251909/193514046488576 

3882036150042265/87747802561511424 

23 

30871937807467/809240558043136 

256340050245714583/6054598376744288256 

24 

123369796036139/3377699720527872 

1538798075750480935/37907050706572935168 
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VII. APPENDIX B: PROOF OF EQS.Qj^ FOR SOME CLASSES OF GRAPHS 


Let 


A(a, V, k) = A^ln(V(j)i!((n — i)!)“) 

with n = [|] and k >2. The bound Eq. Q is A(0, v, k) < 0, the bound Eq. @is A(l,u,fc)<0. 
In the case of polygons we have 


iV(*) = 


V V — i 


V — l \ I 


and thus 


A(a, v,k) = A'^ ^ ( ln(l — —) + ln(l — — ln(l — ^ ^ ) — aln(l — —) 

V V V n 


For a = 0 and all u or o = 1 and v even, 


1 


A(a, V, k) = V — ((i + ly - {2i + IV - (1 - o)2^) < 0 

jr)3 ^ ' 

h^l 


In the case of the complete graphs Ar„, we have 


A(*) = 


v\ 




and thus 


A = A'=-^ ( ln(I - -) + ln(I - - aln(l - -) 

' V V n 


For a = 0 and all 7 ; or a = 1 and even v one gets 


A = A'^-i ^ ^ (-(I - a){2iy - (2* + I)^) < 0 

j=i 


In the case of the complete bipartite graphs Kn,n, we have 


Ni^) ={^]^\ 


and thus 

A{a, V, k) = (2 — a)A^“^ln(I — —) < 0 

n 

Friedland, Krop and MarkstrdmjlQj have computed an average distribution for N{i) in the case of random regular 
bipartite graphs 


Ni^) = 


— i)\ 

(rn)! 


(49) 


and thus 


^ = A‘-‘((2 - a)l„(l - 1) - 1„(1 - ±)) = A‘-‘ y: - (-(2 - a)(ir + j < 0 


J=1 ■ 


The easily derived bipartite mean-field approximation 


/ \ 2 



/r\ 

-1 

. W! 

- 

\^J 

Vn/ 


satisfies the bounds; the proof is similar to that for Kn^n- 
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VIII. APPENDIX C 


Consider the bounds 


< 0 

valid for some fc > 1 and for i > io- Summing from iq to z — 1 one gets 

< A'=-ig(zo) 

Summing again and using 

k\ /n + 1' 


(50) 

(51) 


E 


m) \m + 1 


we can prove by induction that, for h < k 


h-l 


In the case h = k Eq. (52 ) gives 




k-l 


A'=-Vz)<^r A^-'‘+*g(zo) 


(52) 


9{i) < E (* /° )^‘5(*o) 

t=o ^ 


(53) 


for z > zo- 

If a graph satisfies the bound Eq. Q for a given k, then Eq. (53) with g{i) = InP(z), where P is defined in Eq. 


(25), gives an upper bound for N{i). 
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